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Abstract 

In the general matter composition where the multiple scalar fields and the multiple 
perfect fiuids coexist, in the leading order of the gradient expansion, we construct all the 
solutions of the nonlinear evolutions of the locally homogeneous universe. From the mo- 
mentum constraint, we derive the constraints which the solution constants of the locally 
homogeneous universe must satisfy. We construct the gauge invariant perturbation vari- 
ables in the arbitrarily higher order nonlinear cosmological perturbation theory around the 
spatially fiat FRW universe. We construct the nonlinear LWL formula representing the 
long wavelength limit of the evolution of the nonlinear gauge invariant perturbation vari- 
ables in terms of perturbations of the evolutions of the locally homogeneous universe. By 
using the LWL formula, we investigate the evolution of nonlinear cosmological perturba- 
tions in the universe dominated by the multiple slow rolling scalar fields with an arbitrary 
potential. The r function and the potential introduced in this paper make it possible to 
write the evolution of the multiple slow rolling scalar fields with an arbitrary interaction 
potential and the arbitrarily higher order nonlinear Bardeen parameter at the end of the 
slow rolling phase analytically. It is shown that the nonlinear parameters such as Jml, Qnl 
are suppressed by the slow rolling expansion parameters. 
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1 Introduction and summary 



In the inflationary universe scenario, the quantum fluctuations of the scalar fields called 
inflatons are thought to be the origin of the cosmological large scale structures such as 
galaxies and clusters of galaxies. In the slow rolling phase, the quantum fluctuations are 
stretched into the superhorizon scales and stay outside the horizon until they return into 
the horizon in the radiation dominant universe. Therefore the method for investigating 
the evolutions of the long wavelength cosmological perturbations became necessary and 
the LWL method was developed [27| ] [13]] [21]. In the LWL method, we use the LWL 
formulae representing the long wavelength limit of the evolutions of the cosmological per- 
turbations in terms of the quantities of the corresponding exactly homogeneous universe. 
As for the adiabatic modes the exact LWL formula had already been constructed and it 
was used in order to investigate the slow rolling phase [22] ] and the oscillatory phase [12] 
] [6]]. Nambu and Taruya suggested in the multiple scalar fields system also, the LWL 
formula exists [23 ] . Soon in this case the exact LWL formula was constructed [13] ] [21] ] 
and it was used to investigate the multiple oscillatory scalar fields [3] P]. In addition, 
in the paper [13] ] , we presented the flexible way for constructing the LWL formulae in 
the general matter composition. It will be called the Kodama and Hamazaki (KH) con- 
struction in this paper. In the KH construction, the perturbation variables related with 
the exactly homogeneous universe, such as the scalar field perturbation and the energy 
density perturbation, are expressed in terms of the exactly homogeneous perturbations, 
that is the derivative of the exactly homogeneous quantity with respect to the solution 
constant. By solving the spatial components of the Einstein equations, the perturbation 
variables not related with the exactly homogeneous quantities such as vector quantities, 
for example velocity perturbation variables, are expressed in the form of the integral of the 
perturbation variables related with the exactly homogeneous universe which have already 
been determined. We pointed out that the perturbation solution constants contained in 
the expressions of the perturbation variables determined in the processes explained above 
must satisfy the constraint coming from the momentum constraint. The KH construction 
can be applied to the system containing the perfect fluid components having vector de- 
grees of freedom such as the velocity perturbation variables. Therefore by using the KH 
construction, in the most general matter composition where multiple scalar fields and mul- 
tiple perfect fluids coexist, the LWL formula was constructed and used to investigate the 
multiple component reheating and the multiple component curvaton decay [3]. Afterward 
in case of the nonlinear perturbation also, the existence of the LWL formula was suggested 
[15] ] [23 ] . In this paper we give the definition of an arbitrarily higher order nonlinear 
gauge invariant perturbation variables and the exact nonlinear LWL formula representing 
them in terms of the derivatives of the quantities of the locally homogeneous universe with 
respect to the solution constants. 

In order to calculate the present density perturbations and the cosmic microwave back- 
ground anisotropics from the initial seed perturbations, we need to calculate the evolution 
of the scalar fields perturbations on superhorizon scales during the slow rolling phase. In 
the papers [5]] [1], by decomposing the multiple scalar fields into the adiabatic field and the 
entropy fields instant by instant, the evolutionary behaviors of cosmological perturbations 
during the slow rolling phase were discussed. But this study is the local investigation, that 
is, it is based on the Taylor expansion of the evolution equations around the first horizon 
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crossing and the evolutions of the scalar fields in the long time interval were not solved. 
Then in this paper we present the method which makes it possible to trace the evolutions 
of the multiple slow rolling scalar fields with an arbitrary interaction potential for the long 
time period as in the whole slow rolling phase. 

This paper is organized as follows. The section 2 is devoted to the nonlinear LWL 
formula. In the subsection 2.1, under the assumption that the universe is the spatial flat 
FRW universe in the background level, in the leading order of the gradient expansion, 
we give the evolution equations of the locally homogeneous universe. We point out that 
these locally homogeneous evolution equations are similar to the corresponding exactly 
homogeneous evolution equations and that deviations between the two are induced by 
the unimodular factor of the spatial metric 7,^ which has contribution from the adiabatic 
decaying mode. Following the philosophy of the KH construction [I3]], we construct all 
the solutions of the evolution equations of the locally homogeneous universe. We give the 
constraint which comes from the momentum constraint and which must be satisfied by the 
solution constants contained in the solution of the locally homogeneous universe. In the 
subsection 2.2, we give the definition of the gauge invariant perturbation variables of the 
arbitrarily higher order nonlinear cosmological perturbation theory, and prove their gauge 
invariance. We present the nonlinear LWL formula representing the long wavelength limit of 
the evolutions of the gauge invariant perturbation variables in the nonlinear perturbation 
theory in terms of the derivatives with respect to the solution constant of the locally 
homogeneous solutions. The section 3 is devoted to the analysis of the nonlinear evolution 
of the multiple slow rolling scalar fields as an application of the nonlinear LWL formula 
constructed in the previous section. In the subsection 3.1, it is shown that in the slow rolling 
phase the scalar fields evolution equations are simplified by truncation. By estimating the 
truncation error, we establish the accuracy of the truncated evolution equations in the 
slow rolling expansion scheme. In the subsection 3.2, the ideas of the r function and 
the potential are introduced and they are shown to enable us to trace the multiple 
slow rolling scalar fields in the whole slow rolling phase analytically. By adopting the 
T function as the evolution parameter, the truncated evolution equations of the multiple 
slow rolling scalar fields are simplified enough for their solutions to be written analytically. 
From the scalar fields solutions we can easily calculate the N potential which allows us 
to calculate the arbitrarily higher order Bardeen parameter at the end of the slow rolling 
phase from the initial scalar fields perturbations at the first horizon crossing. In the 
subsection 3.3, by the r function and the potential we calculate the various perturbation 
variables such as the Bardeen parameter, the entropy perturbations, the gravitational wave 
perturbation and their spectrum indices in the case of the multiple quadratic potential 
whose truncated evolution can be exactly solved. In the subsection 3.4, we consider the 
effect of the interaction between multiple slow rolling scalar fields in the case where the 
masses of the scalar fields do not satisfy any resonant relations. We point out that this 
problem is related with the well known Poincare theorem about the linearization. In 
the concrete model, we calculate the various quantities such as the potential and the 
nonlinear parameters /nl, Qnl introduced in the paper [H]]. In the subsection 3.5, we 
investigate the resonant interaction between the multiple slow rolling scalar fields. We 
show that the potential has no singular part by introducing the resonant interactions. 
The section 4 is devoted to the discussions. The appendices are devoted to the proofs of 
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the propositions presented in the main content of this paper. 

§2 Nonlinear LWL formula 

In this section, we derive the nonhnear LWL formula in the most general matter composi- 
tion whose energy momentum tensor is divided into A = {S, f) parts where S represents 
Ns components scalar fields (pa {a = 1, 2, ■ ■ - , A^^) and / represents Nf components perfect 
fluids pa (a = 1,2,- ■ -jNf). The content of this section is the nonlinear generalization 
of the content of the paper P]. Our results in this section can be applied not only to 
the multiple slow rolling scalar fields, but also to the multicomponents reheating and the 
multicomponents curvaton decay [20|] P]. 

In the subsection 2.1, we give the evolution equations of the locally homogeneous uni- 
verse in the leading order of the gradient expansion and we construct all the solutions of 
these locally homogeneous universe evolution equations. In the subsection 2.2, we give the 
definition of the nonlinear gauge invariant perturbation variables and we give the nonlinear 
LWL formula representing the long wavelength limit of the evolutions of the gauge invari- 
ant perturbation variables in terms of the derivative with respect to solution constant of 
the evolutions of the locally homogeneous universe determined in the subsection 2.1. 

Our theory has two small expansion parameters: One is e characterizing the small 
spatial derivative, that is, the small wavenumber, and the other is 6c characterizing the 
amplitudes of the higher order perturbations. In general, under our present scheme, all 
the terms are classified as 0(e^(5') where k, I are appropriate nonnegative integers. Since 
we are interested in the full nonlinear perturbations, in the subsection 2.1 we will not 
Taylor expand with respect to 6c- But since we treat the evolutions of the long wavelength 
perturbations only, we will expand with respect to e and we will drop terms which are 
small compared to the leading order by O(e^) order quantity. This process corresponds 
to the leading order of the gradient expansion [22] ] [ISl ] [2S] ] and the universe treated in 
this way is called the locally homogeneous universe. The locally homogeneous evolution 
equations are very similar to the exactly homogeneous evolution equations. The evolutions 
of locally homogeneous physical quantities which have counterparts in the corresponding 
exactly homogeneous universe can be determined as easily as the exactly homogeneous 
evolutions are determined. This is the attractive point of our LWL method. 

In our scheme, we can assign dt^ij = 0{6c), d^^ij = 0{6c) where is the unimodular 
factor of the spatial metric defined by (12. 5p . since we require only that the universe should 
be the spatially fiat FRW universe in the background level. So in the leading order of 
the gradient (e) expansion the terms containing dtjij, df'jij cannot be dropped. But in 
the papers [15] ] [26] ] , mainly in order to avoid the computational complexity, the authors 
assigned dt^ij = O(e^), d^^ij = O(e^) and discarded the terms containing dtjij, d^jij. 
But this is confusion between the different small parameters e, 6c and cannot be justified. 
In this section, under the more natural assumption 'jij = 6ij + 0{6c), we show that the 
evolution of jij can be solved analytically and that the evolution equations of the other 
dynamical variables are simple enough to be solved analytically in spite of inclusion of 
the terms containing dt'jij, df^ij. By doing so, in the leading order of the gradient (e) 
expansion, the evolutions of all the modes 2Ns + 4:Nf + 4 are obtained consistently in the 
universe where Ns scalar fields and Nf perfect fluids coexists. In the papers [15]] [26]], 
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the evolutions of the several dynamical variables become higher order O(e^), therefore the 
initial conditions of such variables are unnaturally constrained to too small quantities. 



2.1 Evolution of the locally homogeneous universe 

First we present the Einstein equations G^^i, = n'^T^u where = SuG is the gravitational 
constant, based on the 3 + 1 decomposition based on the paper |25|]. The metric 

ds^ = g^^dx'^dx'', (2.1) 

is written by 

goo = + (2.2) 

goi = A, (2.3) 

9ij = 7*i, (2.4) 

where a is the lapse and f3i is the shift vector and /3* := Y'^Pj- Greek indices take the 
values /i, z/ = 0, 1, 2, 3 and Latin indices take values i,j = 1, 2, 3. The spatial metric •jij is 
decomposed as 

7ij = a'^lij, det(%) = 1, (2.5) 

where a can be interpreted as the nonlinear generalization of the scale factor. The extrinsic 
curvature of the t = const hypersurface is given by 

= ^ (7., - A/?, - I^iA) , (2.6) 

where Di is the covariant derivative with respect to •jij. The extrinsic curvature is decom- 
posed as 

= ^^ijK + a^Aij, (2.7) 
Y'Aj = 0. (2.8) 



Then we obtain 



a \ a 

The energy momentum tensor is given by 



-K=-[3-- ) . (2.9) 



T^u = ip + PWu, + Pg^,, (2.10) 

where p, P, are the energy density, the pressure, and the 4 velocity of the total system, 
respectively. The 4 velocity is written as 



■° - [a' - iP, + V,) + v')]-'^\ (2.11) 
u' = u%\ (2.12) 
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where Vi is the 3 velocity of the total system and := yhjj. By using := {—a, 0, 0, 0) 
which is the unit vector normal to the time slices, the 3 + 1 decomposition of the energy 
momentum tensor is given by 



E 
J. 



The Hamiltonian and momentum constraints are written as 



DiA\ - -D,K 



(2.13) 
(2.14) 

(2.15) 



(2.16) 
(2.17) 



where the indices of A^j is raised by 7*-^ which is the inverse matrix of 7^. The evolution 
equations for are written as 



{^^-|3''^k)'^^J = -2aA,+7,fca,/5'^ + 7,fca#-^7^A/5'■ 
The evolution equations for K^j are given by 



-D^D^a + ^a{E + S\) , 



(2.18) 
(2.19) 



(2.20) 



a I R 



{dt-(5^dk)A,^ = - 

Qj 

+a (KAi^ - 2AfeA^.) + A,kdj(5'' + A^A(5 



2 T 



(2.21) 



Rij is the Ricci tensor of the metric and R = Y^Rij, S\ = 'y'^''Ski- 

Next under the 3 + 1 decomposition, we rewrite the evolution equations of the energy 
momentum tensor T^^^ of the A component: 



(2.22) 



where A := {a, a), a is the scalar fields index and a is the perfect fluids index. The energy 
momentum transfer vector Qa// is decomposed as 



Qah — QaUh + fAfj., 

u''fA^. = 0, 



(2.23) 
(2.24) 
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where Qa, Jau are the energy transfer, the momentum transfer vector of the A component, 
respectively. The energy momentum tensor of the total system can be expressed as the 
sum of the each component energy momentum tensor T'^^: 

(2-25) 

A 

The energy momentum conservation, that is, the Bianchi identity gives 

5^g^^ = 0. (2.26) 

A 

We consider the scalar field component = {(pa)- The energy momentum tensor of the 
scalar fields part is given by 

(r(;)5 = V^c/) • - ^ [g'"'V,<l> ■ + 2U\ g,, (2.27) 

The energy momentum tensor of the scalar fields part {T^i,)s can be written as the perfect 
fiuid form by identifying 

PS = —g'^Vpck-VA + U, (2.28) 



= -^/"V,</)-V,0-t/, (2.29) 



sign(0„) {-gp''Vp<t> ■ Va<t>) ' 

where the minus sign of is adopted by requiring UaQ — —ol + O(e^) in the gradient 
expansion scheme, where e is the small parameter characterizing the spatial derivative 
defined below. When we assume that the energy momentum transfer vector of the scalar 
fields part {Qfj.)s is given by 

(Q^)5 = SaV^cPa, (2.31) 

where the source function Sa describes the energy transfer from the scalar field (pa to other 
components, the evolution equation of the scalar fields components V,/(Tp5 = {Q,j)s holds 
if the scalar field (pa satisfies the phenomenological equations of motion of the scalar field 

(pa- 

dU 

d(pa 

where 

□0a = {-gr''^d,[{-gf'^g^^d,(P^, (2.33) 
g := det(^^,). (2.34) 

Since we are interested in the cosmological perturbations on superhorizon scales, we 
put the gradient expansion assumption defined by 

9, = 0(e), A = 0(e), = 0(e), /a^ = 0(e), (2.35) 
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where e is the small parameter characterizing the small wavenumber of the cosmological 
perturbations. By the assumption di = 0(e), we assume that the spatial scale of all 
the inhomogeneities is of the order of 1/e, that is, all the physical quantities which are 
approximately homogeneous on each horizon can vary on the superhorizon scales. The 
local homogeneity and isotropy in the horizon guarantee that the vector quantities such 
as Pi, Vi, fAi are of the order of 0(e). Unlike the papers [15]] [26|] where dt^ij = O(e^) 
is assumed, from the requirement that our locally homogeneous universe should be the 
spatially flat FRW universe in the background level, we assume that 



lij = Sij + 0(4). 



(2.36) 



where 6c is the small parameter characterizing the higher order perturbations. Since we 
consider the leading order of the gradient (e) expansion without expanding with respect 
to 6c, we must keep the terms containing dt^ij = 0{6c), d^^ij = 0{6c) as explained in the 
beginning of this section. 

Under the gradient (e) expansion scheme, the relations between the total system quan- 
tities and the each component quantities are given by 



P 
P 

h 

hvi 







A 
A 



A 

J2hAVAi + 0{e^ 

A 

J2qa, 

A 

fAi, 



(2.37) 
(2.38) 
(2.39) 
(2.40) 
(2.41) 
(2.42) 



A 



where Ha is the A component enthalpy defined hj Ha '■= Pa + Pa- 
As for the fluid component a, V^T^q = Qao gives 



Po 



-3H{pc, + Pa) + Qaa + 0(e^ 



where H is the Hubble parameter defined by H := a/a, and integrating V^T^j 
respect to t gives 



a-^ a-^ 



dtaa 



to 



-diPa - -Diaha + Qc 
a 



o(a, 



(2.43) 
Qai with 

(2.44) 



where to is the initial time and C^i '■= Cai{x) are the integration constants. 

As for the scalar field components, from (12.321) . the equation of motion of the scalar 
field (pa is given by 



1 



K + S- 



— < 

a 



dU ^ 



(2.45) 
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Since the source function Sa is the scalar quantity, we can assume that Sa is the function 
of other scalar quantities. As such scalar quantities, we can adopt (pa, 

T2a := sgn(9o0a) (-/'^Vp0„V.0a)'^' = -0a + ^(e^), (2.46) 

a 

or 

T2a ■■= U^'^^<Pa = -k + 0{^) (2.47) 

a 

where ti^ is the arbitrary unit time like vector field and 

T3:= VX = 3-- + 0(e2). (2.48) 
a a 

Therefore in the leading order of the e expansion, the form of the source functions Sa can 
be given by 

Sa = Sa (<Pa, -<Pa, S--) + ©(e^). (2.49) 

For example, 

Sa = Ta-k, (2.50) 

a 

where Ta is the decay constant of the scalar field (pa, is the most simple source function. 

Since until now we have presented the leading order of the gradient (e) expansion of 
all the locally homogeneous evolution equations, we will construct all the solutions of the 
evolutions of the locally homogeneous universe. From fl2.19p . we obtain 

By substituting the above equation into f l2.2ip . we obtain 



Mi" 



i^hi^M + MM-^M + Oie^), (2.52) 
\dt ) 

where M := {^ij). By neglecting O(e^) order terms, we obtain the solution as 

M = Rexp 



dt^T 

to « . 



(2.53) 



where R = R{x), T = T{x) are time independent matrices. Since M is a unimodular 
symmetric matrix for an arbitrary t, R is unimodular symmetric, T is traceless, and RT 
is symmetric. By using (12.531) . we obtain 

Ai.A'^ = Artr (mM-^MM-^) = (2.54) 

CT := ^tr(T2), (2.55) 

where ct = ct{x) is a time independent constant. By using the above results, (12.161) gives 
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(K20^ gives 

By eliminating H in (12.561) and (12.571) . we obtain the well known continuity equation of 
the total system in the expanding universe as 

p = -3H{p + P). (2.58) 

Eqs. (12.431) (12.451) and the above three evolution equations agree with the exactly homoge- 
neous evolution equations if the proper time slicing a = 1 [10 ] and ct = 0. But we cannot 
assume ct = since the solution constants must satisfy the constraint originating from the 
momentum constraint (12.17^ as 

= 0. (2.59) 

to 

For the derivation of h2. 59\) . see Appendix A. 

We consider the long wavelength limit of all the solutions of the evolution equations 
of the locally homogeneous universe. The unimodular factor of the spatial metric jij is 
obtained by (12.531) . These induce the deviation between the true locally homogeneous 
universe and the corresponding exactly homogeneous universe, for example ct terms in 
(I2.56P and (12.571) . They contain the scalar adiabatic decaying mode which was carefully 
treated in the papers [T3| ] ^ ], since in the scalar part in the linear cosmological per- 
turbation theory it induces the deviation between the long wavelength limit of the true 
perturbation solutions and the derivative of the exactly homogeneous universe with respect 
to the solution constant. According to the philosophy of the KH construction [I3]], the 
expressions of pa, 4>a related with the exactly homogeneous quantities are obtained by solv- 
ing the evolution equations (12.43^ (I2.45P under the Hamiltonian constraint ( I2.56P and the 
proper time slicing a = 1. Under the proper time slicing a = 1, (I2.43p . (I2.45P and (12.561) 
are almost the same as the counterparts of the exactly homogeneous universe. Therefore 
solving the former true locally homogeneous evolution equations requires as little labor as 
solving the latter exactly homogeneous evolution equations. As for the velocity perturba- 
tions of the fluid components ha{Pi + f«j) which are the vector quantities not related with 
the exactly homogeneous quantities, their evolutions are given by (12.441) in whose right 
hand side the second integral term contains Pa, ha, « which have already been determined 
by the previous process. The solution constants must satisfy the momentum constraint 

Let us count the degrees of freedom. As for (12.530 . R is unimodular symmetric, T is 
traceless, and RT is symmetric, therefore R, T have 5 degrees of freedom, respectively. But 
by the coordinate transformation = /*(x) {i = 1, 2, 3) the 3 degrees of freedom of R can 
be made vanishing. According to (I2.43P (I2.45p . the densities pa and the scalar fields (pa have 
Nf, 2Ns degrees of freedom, respectively. According to (12.440 . the fluids velocities Vai have 
3A^/ degrees of freedom. The momentum constraint (I2.59P gives 3 constraints. Therefore 
the total degrees of freedom is 5 + 5- 3 + A^/ + 2A^5 + 3A^/ - 3 = 4 + 4A^/ + 2Ns. Then 
we have obtained all the solutions of the evolution equations of the locally homogeneous 
universe in the leading order of the gradient (e) expansion. 
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From the time dependence of all the solutions, we can interpret the physical roles of all 
the solutions. 2 from R can be interpreted as the gravitational wave growing modes. 5 from 
T can be interpreted as the 2 gravitational wave decaying modes, the 1 adiabatic scalar 
decaying mode and the 2 adiabatic vector decaying modes. By the 3 momentum constraints 
fl2.59p . 3 of Cja are adjusted. In the remaining (3A^/ — 3) CjaS, the {2Nf — 2) entropic 
vector decaying modes and the (Nf — l) entropic scalar decaying modes are contained. The 
Nf densities pa have the Nf scalar growing modes and the Ns scalar fields (pa have the Ns 
scalar fields growing modes and the Ns scalar fields decaying modes. 

2.2 Gauge invariant variables and the derivation of the LWL for- 
mula 

In this subsection, we give the definitions of the gauge invariant perturbation variables 
in the arbitrary higher order perturbation theory in the leading order of the gradient 
(e) expansion and the LWL formula representing the evolutions of these gauge invariant 
perturbation variables in terms of derivative with respect to the solution constant of the 
corresponding physical quantities of the locally homogeneous universe. 

As for M = {jij), since the homogeneous parts of R, T defined by (12.531) are deter- 
mined by the fact that the background metric is the spatially flat FRW universe, R, T are 
expanded as 

R,j{x) = + (2.60) 

k=l 

oo _ 

^^^•(^) = (2.61) 

k=l 

where 6 in the above is the operator generating the higher order perturbation quantities. 

As shown in the previous subsection, in the leading order of the gradient (e) expansion, 
the evolution equations of the locally homogeneous universe which have counterparts in 
the evolution equations of the exactly homogeneous universe do not contain the spatial 
derivative. Therefore in the expression of the solution of the evolution of an arbitrary 
locally homogeneous physical quantity (related with the exactly homogeneous quantity) A 
such as the scalar fields (pa and the fluid energy densities pa, all the dependences on the 
spatial coordinate x are contained in the (time independent) spatial dependent integration 
constants C{x): 

A = A{t,Cix)). (2.62) 

Since the solutions of the evolution equations of the locally homogeneous universe contain 
the nonlinear effect in the full order, the locally homogeneous physical quantity A can be 
expanded as 

oo _ 

A(t,a;) = A(t) + ^-5'^A(t,x). (2.63) 

k=l 

The higher order perturbation effects are induced by the dependences on the spatial coor- 
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dinate x. Therefore the each solution constant Ci{x) can be expanded as 



a{x)=C, + Y,y''C,{x), (2.64) 

k=l 

whose background part is spatially independent. The fc-th order nonlinear perturbation 
5^ A can be expressed in terms of the nonlinear perturbations of the spatial dependent 
integration constants 5^C{x). In this paper, the expressions representing 6^ A in terms of 
5^C{x) are called the LWL formula. In order to derive the LWL formula, we propose a 
simple mathematical trick. In the leading order of the gradient (e) expansion, the expression 
of the solution of an arbitrary locally homogeneous physical quantity A can be written as 
function of the time t and the integration constants C as shown in fl2.62p . We assume 
that all the integration constants C depend on one parameter A imaginarily instead of x. 
As understood below, A represents the spatial coordinate dependence symbolically. The 
physical quantity A can be expanded as 



°° 1 



A:=0 



(2.65) 

A=0 



The full order nonlinear solution is formally recovered by setting A = 1. Since both the 
perturbation 5 and the A differentiation d/d\ are the derivative operators satisfying the 
same chain and product rules, and the fc-th order A differentiation d^ / c/A'^ - ■■ | a=o is multiplied 
by A'^, we can use the A differentiation to track the algebraic behaviors of the perturbation 
5. By comparing (12.651) with (12.631) . we can read the correspondences given by 

A{X) 



d\^ 



^S^Ait.x), (2.66) 

A=0 

The gauge transformation can be expressed by the Lie derivative L{T) as 

A(A, /i) = exp{/iL(T)}A(A, = 0), (2.67) 

where /i is the parameter characterizing the size of the gauge transformation and v4(A, ^ = 
1) is the transformed variable and A{X, n = 0) = A{X) is the original variable. When the 
vector field T := T'^d^ in the Lie derivative L{T) can be expanded in terms of A, the zeroth 
order term T(A = 0) is zero since we consider the infinitesimal gauge transformation. By 
differentiating (12.671) with respect to A, afterward putting A = 0, we obtain the well known 
expressions of the gauge transformations of 6'' A [3]]: 

6~A = L1A + 6A, (2.68) 
6^A = L2A + LiLiA + 2Li5A + 5^A, (2.69) 

where 6'^A is the gauge transformed perturbation variables of 6"'A, and Lk is the Lie 
derivative induced by S'^T; Lk := L{5'^T). The gauge transformation law (I2.67P is the 
solution of the differential equation: 

|-A(A,/i) = L(T)A(A,/x), (2.70) 
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which is much simpler than the individual gauge transformation expressions of S^A. As 
for the vector field T := T^^d^ inducing the Lie derivative L{T), we can assume that 

T = 0(e), (2.71) 

since we consider only the gauge transformations keeping the local homogeneity and isotropy 
in the horizon. 

By using the differential equation (12.701) . we can show the following proposition: 

Proposition 1 For an arbitrary scalar quantity A, the perturbation quantities D^A, 
D'^{A/a) where D is defined by 

d\ d\ a dt 

are gauge invariant up to order O(e^) error. For the lapse function a, An defined by 

A d 1 da\^ (2 73) 

" \dX dt ddXj 

is gauge invariant up to order 0{e^) error. 
For the proof, see Appendix C . 

D"'A, D"'{A/a) and An are higher order generalization of DA, DA and A defined in 
the paper P], respectively. D"-{A/a) and An are not independent. For example 




^{(DA)-a-M}. (2.74) 



Therefore by putting A = (j)a, A = pa, we can use D^'cpa, -D"(0a/«), D"-pa as the indepen- 
dent perturbation variables. Our DA is almost the same as the well known gauge invariant 
perturbation variable, for example, for A = (pa the Sasaki Mukhanov variable in the linear 
perturbation theory [10]] [H]] [IS] P]. Our D'^A agrees almost perfectly with the gauge 
invariant quantities introduced by the second order gauge invariant perturbation theory 
[16] ] [T7] ] [19]]. The reason of the tiny deviation between our gauge invariant perturbation 
quantity and the gauge invariant perturbation quantity of the first and the second order 
gauge invariant perturbation theory is that we assume T* = 0(e) based on the gradient 
expansion scheme. In our scheme fl2.7ip . we can regard the scale factor a as the scalar 
quantity up to O(e^). Therefore in the same way as in the proof of Proposition 1 we can 
show that 

Cn := 5" In a, (2.75) 

d\ d\ p dt ^ ^ 

are gauge invariant up to O(e^). (n is the higher order generalization of the well known 
Bardeen parameter [2]] [28]]. 

We discuss the influence of the replacement of the evolution parameter. Concretely 
we consider replacing the old evolution parameter t with the new evolution parameter as 
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the scale factor a. With the scale factor a as the evolution parameter, the D operation 
defined by fl2.72l) can be written in the more simple form. d/dX in fl2.72p is the A derivative 
with if: fixed, that is, operating on the locally homogeneous physical quantity such as Pa, (pa 
expressed by using t as the evolution parameter. However an arbitrary locally homogeneous 
quantities can be also expressed by using the scale factor a as the evolution parameter 
instead of t. In this case, we can consider the differentiation {d/dX)a taken at fixed a. The 
differential operators in the two group {d/d\,d/dt) {{d/dX)a, {d/da)a) are commutative in 
the each group, but the differential operators belonging to the different groups, for example 
d/d\ and {d/da)a, are not commutative. By using the relation 

4^=(r]+^(l), (2.77) 



dX \dX J „ dX \da 
the D operation f l2.72p can be expressed much more simply as 



D=[4t] ■ (2-78) 



dX 

Therefore D"'A, D"'{A/a) can be written in the very simple way as 

D''A=(^ \ A, D"{-\ = [—-] {-]. (2.79) 



dX 



a 





In this way, in the LWL formalism where the scale factor a, not t, is used as the evolution 
parameter |7|] [8]] [9|], the expressions of the solutions of D'^A, D"'{A/a) can be obtained 
by calculating only a single term written with the higher order {d/dX)a derivative of the 
solution of the corresponding locally homogeneous physical quantity. In the same way, 6 
defined by (12.761) can be expressed in the word of the scale factor a as 



dX J ^ \\dX J \da J „ \ \da 

\ /a v\ /a \ / a y \ / u, 

As for the Bardeen parameter defined by f l2.75p . we can show the following propositions. 

Proposition 2 Ij is conserved for arbitrary values of integration constants C{X = 0), 
all (n (n > 2) are also conserved, and they can be expressed as 

d " 



Cn = ( ^ ) Cl (2.81) 



a 



For the proof, see Appendix D 

Proposition 3 When P = P{p) holds, (n (n > 1) can be expressed as 

1 



1 / 9" 

Cr 



3 \dX 

and all (n (n > 1) are conserved. 



dp 



P + P{p) 



(2.82) 
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For the proof, see Appendix E 

All the definitions of the nonlinear gauge invariant perturbation variables in this sub- 
section are written in terms of the A differentiations. As explained at the beginning of 
this subsection, the A differentiation implies not only the symbol of the higher order non- 
linear perturbation of the physical quantity, but also the process of taking the derivative 
with respect to the integration constants of the corresponding locally homogeneous quan- 
tity. Therefore all the definitions of the nonlinear gauge invariant perturbation variables 
in terms of the A differentiations can also be regarded as the LWL formulae themselves, so 
from now on they will be called the LWL formulae. 



§3 Nonlinear evolution of the multiple slow rolling 
scalar fields 

In this section, as the application of the LWL formula derived in the previous section, 
we consider the evolutions of the long wavelength nonlinear cosmological perturbations 
in the universe dominated by the multiple slow rolling scalar fields. The r function and 
the potential introduced in this section are the useful tools for tracing analytically the 
evolutions of the multiple slow rolling scalar fields in the long time interval. We calculate 
spectral indices of the linear cosmological perturbations. In the interacting system, we 
derive the formulae giving the amplitudes of the nonlinear Bardeen parameters at the end 
of the slow rolling phase in terms of the initial scalar fields perturbations, and calculate 
the nonlinear parameters /tvl, Qnl P3!] representing the non-Gaussianity of the Bardeen 
parameter. 



3.1 Evolution of the multiple slow rolling scalar fields 

In the slow rolling phase, the scalar fields 0^ roll slowly on the potential U . The potential 
energy U which hardly changes triggers the exponential expansion of the universe. The 
Hubble parameter is large compared to the masses of the scalar fields. The ratio of the 
kinetic energy part in the whole energy density p is small compared to the contribution 
from the potential energy U . In the investigations of the evolutions of the scalar fields under 
this situation, it is effective to use the transformation by which the evolution equations of 
the slow rolling scalar fields system are greatly simplified, that is the effects of the time 
derivatives of the scalar fields pa '■= (pa/oi on the evolutions of the scalar fields (pa are 
eliminated. 

In this section, we consider the multiple slow rolling scalar fields (pa under the conditions: 



= 0, U 



+ (3.1) 



where U\nt is the sum of m-th order monomials (m > 3). As the independent variables, we 
adopt 0a and pa (pa/ot- By nondimensionalizing the dynamical variables as 

a (pa Pa Ct 

^ a, ^ (pa, >Pa, ^ Ct, (3.2) 

ao 00 Po cto 
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and the parameter as 

— rria, (3.3) 
mo 



that is 

P 



we obtain the dimensionless parameters: 

e.:=—. V':=^, ■.= -^^-. (3.5) 

e* is the small constant representing the ratio of the mass scale to the Hubble parameter. 
This e* is different from e characterizing the small wavenumber in the previous section. 
rf is positive constant since we consider only the scalar fields with positve definite kinetic 
parts. We assume that e,,, ~ ?7 ~ z/ <^ 1. Then the evolution equations of 0a, Pa and ct are 
given by 

Ct = const, (3.8) 

where the evolution parameter is the e- folding number of the scale factor := In a and 
the energy density p is given by 



P 



yT^p' + U + u'ct^. (3.9) 



By using p^a^ defined by 



(1) , 1 1 dU 



P^a^ ■■=Pa+--^J^, (3.10) 

which represents the deviation of Pa from the truncated slow rolling solution, the evolution 
equations fl3.6pfl3.7p can be rewritten as 

= v'Fa{<l>) + r]'U<P,p^'\cT,N), (3.11) 

J^P^a^ = ~3p^^^ + v'9a{<P,P^'\cT,N) (3.12) 

where 

-Ht^It' (3.13) 
3 r/^ (/ 0(pa 

and 

\fa\<\p^'^+V^ \9a\<l, (3.14) 

for an appropriate complex domain containing the real interval where we consider the 
motion of 0^, p^\ N. In this section, in all the inequalities we omit all the finite constants 
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and \p^^^\ is interpreted as the quantity bounded by M|p*^^^| for some positive constant M 
and 

\P^'^\--=J2\Pa^\- (3.15) 

a 

From now on, we will simply write p^^ as pa for notational simplicity. In this section, we 
consider the evolution for < < l/r/^ during which the scalar fields roll slowly on the 
potential U. For a function f{N), let us define ||/|| by 

(3.16) 



= sup \f{N)\. 

0<Af<l/»?2 

Under these notations, following propositions hold. 

Proposition 4 Let k nonnegative integer, and 6c small positive constant. Under the 
initial conditions 



Qk 



0(0), 



Qk 



^P(O), 



Qk 



for < N < 1/t]'^, the upper bounds of the independent variables are given by 



Qk 
Qk 



p 



< 5l 



-3N 



(3.17) 

(3.18) 
(3.19) 



For the proof, see Appendix F. 

In the above and from now on, as for A differentiations of the physical quantities at 
the initial time N = 0, the suffixes a implying "a fixed" are omitted since what the A 
differentiations operate on do not contain any a = dependent parts. 



Proposition 5 Let k nonnegative integer. The differences A(f)a '■= <pa — (t>a where (j) obey 
the exact evolution equations ^3.11\) ^3.12\l and cf) obey the truncated evolution equations as 



are bounded as 



Qk 



A(j) 



Qk 



under the initial conditions 



Qk 



A(/)(0) = 0. 



(3.20) 



(3.21) 



(3.22) 



For the proof, see Appendix G. 

According to Proposition 5, if we want to investigate the evolution of the leading order 
in the slow rolling (77^) expansion, we only have to solve rather simple evolution equations 
fl3.20p . In the following subsections, we study the evolutionary behaviors of the above 
evolution equations (13.201) . 
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3.2 the r function and the potential 

In this subsection, we introduce the r function and the potential which enable us to 
trace the multiple slow rolling scalar fields in the long time interval analytically. By using 
the r function as the evolution parameter, the evolutions of the scalar fields (f)a and the old 
evolution parameter N can be expressed in the form of the simple analytic function of r. 
From these expressions of (pa in terms of the r function, we can calculate the potential 
in the simple way. The potential is written in terms of the initial values of the scalar 
fields 0a(O) only. It not only represents the difference between the e-folding number of 
the scale factor at the end of the slow rolling phase and that at the first horizon crossing, 
but also has the complete information as to the nonlinear curvature perturbations. The A 
differentiations of the N potential generate the S formulae connecting the amplitudes of all 
the higher order Bardeen parameters at the end of the slow rolling phase with the scalar 
fields perturbations at the first horizon crossing. 

Under the conditions (13.11) . we will investigate the evolutionary behaviors of the solution 
of the evolution equations which are obtained by the truncation explained in the previous 
subsection. 

From the present subsections, we will call off the nondimensionalization in the previous 
subsection, since the truncated evolution equations fl3.23l) have already been made suffi- 
ciently simple. In the multiple scalar fields case, the evolution equations (13.231) cannot 
be solved analytically, and the scalar fields 0a cannot be expressed in the form of the 
well known function of A^. So by replacing the old evolution parameter A^ with the new 
evolution parameter r, we decompose (I3.23P into two parts: 

- -w: <-^' 

= k'U. (3.25) 

ar 

The new evolution parameter r introduced in the above equations will be called the r 
function from now on. By introducing the r function, the evolution equations become 
simple enough to be solved analytically. In the potential (13.11) . we can easily solve (13.241) 
by iteration, then we can get 0a (t) expressed in the form of the analytic functions of the r 
function. By substituting these (pai.'^) into U in (I3.25p . and integrating (I3.25P with respect 
to the r function, we obtain the expression of the old evolution parameter A^ in terms of 
the new evolution parameter r: 

N= ! drn'^U. (3.26) 



Jo 

The expressions of (pa, N in terms of the r function describe the dynamical evolutions of 
our multiple slow rolling scalar fields system completely. 

From this subsection, we adopt the r function as the evolution parameter. We introduce 
{d/dX)r as the A differentiation taken at the fixed r, that is, operating on the locally 
homogeneous quantities expressed by using the r function as the evolution parameter. In 
order to exaggerate the fact that the r derivative and {d/d\)r are commutative, we use 
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{d/dT)r as the r derivative from now on. By using {d/d\)r, ((9/5r)^, D defined by (12.721) 
can be expressed as 

0-(l)-im (3.27) 



dX J \\dX J I Kdr J \ Kdr 
and 6 defined by (12.7611 can be expressed as 



(3.28) 



dX J \\dXj I \dr ) \ \dr 

wfiere we used 

d f d \ dr f d 



dx \dxj^^ dx\dT)^' ^^-^^^ 

and p = U which holds under the present truncation (I3.23p . Then by using the r function 
as the evolution parameter, the Bardeen parameter C„ defined by (I2.75P can be decomposed 

as 

Cn = ^iV--5M(0,0). (3.30) 

By N, we represent 

/■oo 

N ■= dTK% (3.31) 
Jo 

and this N will be called the potential from now on. In the above and from now on, 
as for A differentiations of the physical quantities at the initial time r = 0, the suffixes r 
implying "r fixed" are omitted since what the A differentiations operate on do not contain 
any r dependent parts. A{2n, k) is defined by 

oo 



A(0,0) := 4 / drU, (3.32) 

ly f&^d^ 

2/ V 5A^ 



A{2n,k) := (-- U-^l AM. (3.33) 



By using A{2n, k), 6 can be expressed as 



^-il)^l^(l)- (3.34) 



dXj^ 2A(4,0)V5r^ 

Afterward we will prove the fact that as for the Bardeen parameter (n at the end of 
the slow rolling phase the first term in (13.301) has leading contribution in the slow rolling 
expansion scheme. In order to prove this statement, we put forth several assumptions. All 
the investigations in this section will be established under the following assumptions: 

(i) All the masses of the scalar fields are of the same order. 

ml ~ m^, (3.35) 

where 

:= min{m^}. (3.36) 
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(ii) The interaction potential Ui^t is the sum of the m-th order monomials (m > 3) and 
satisfies 

\Uint\<ficm^^\M', (3.37) 
for 101 < 00- 00 is defined by the positive constant of the order of 

00 ~ 0a(O), (3.38) 

and 101 is defined by 

|0p:= J]|0.r, (3.39) 

a 

and /ic is the small positive constant characterizing the interaction strength. 

(iii) The n-th order perturbation at the initial time r = is of the order of 

Qk 

Q^MO) ~ ^,'00, (3.40) 

where 6c is the small constant characterizing the perturbation size. 
Under these assumptions, the following proposition holds. 

Proposition 6 The following estimations hold. 

= tS^['/>^(0)]+<-Vc«>, (3.41) 

a 

A{2n, k) = J2 HT ^[<^a(0)] exp [-2m^]+ < /x.m^-^^V^ exp [-2mV] >,(3.42) 

a 

where < M > is the quantity bounded by M. 
For the proof, see Appendix H . 

At the end of the slow rolling phase when the Hubble parameter is of the order of 
the scalar fields mass, the values of the scalar fields are of the order of the Planck mass 
0a(T) ~ l//t. According to Proposition 6, at the end of the slow rolling phase, we obtain 
the estimations as 

^'iV ~ (3.43) 



, 2 

,2n(:fe ,2 ' 



A{2n,k) ~ , (3.44) 

where Hq is the Hubble parameter at the initial time r = 0. By using the above estimations, 
the second term of (13.301) can be estimated as 

4?M,0,0)^K^*>?(^)^|LiV.(|^)\ (3.45, 
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Since the second term of (13.301) is suppressed by the slow roUing parameter [m/HQ)'^ com- 
pared to the first term, the main part of the nonhnear n-th order Bardeen parameter at 
the end of the slow rolling phase is given by the A derivative of the potential: 

(n = g^N. (3.46) 

So in order to obtain the final amplitude of the Bardeen parameter we have only to 
calculate the potential N. 

In particular case, the potential can be calculated very easily. When the scalar fields 
potential is written in the separable form 

U = J2UaM, (3.47) 

a 

the potential can also be expressed in the separable form as 

N = Y,Nia), (3.48) 

a 

where 

rMo) / /BIT \ 
N{a):=K^ rf0 Jf/ . (3.49) 



3.3 Exactly solvable model; noninteracting case 

In this subsection, in the multiple free fields case, we give the Bardeen parameter and 
the entropic perturbations and the gravitational wave perturbations and their spectral 
indices. Since we can obtain the analytic expression of 0a (t), all the results obtained 
in this subsection have no black boxes originating from the S formulae connecting the 
amplitudes of the final physical quantities and those of the initial physical quantities, 
which were unknown in the paper [H ] . 

We consider the exactly solvable model defined by 

a 

and calculate the various physical quantities. The A^ potential is calculated as 

^=tS<^'(0)' (3.51) 

a 

then the Bardeen parameter C„ at the end of the slow rolling phase is calculated by fl3.46p . 
From now on, we consider the linear perturbations, and calculate their spectral indices. 
We consider the entropy perturbation between (pa and 4^, defined by 

Sa, ■= -3H (-Dpa - -Dp,] . (3.52) 

\Pa Pb J 
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In the present case (13.501) . Sab is given by 

In order to calculate spectral indices, we have to know the amplitudes of the quantum 
fluctuations of the scalar fields and the gravitational waves at the first horizon crossing 
r = 0: 



dX dX 



K^W dX dX 



H'6ab ~ K'U6ab, (3.54) 

H'^ ~ K^U. (3.55) 



When we calculate the above correlation functions, for simplicity we do not take into 
account the first order slow rolling corrections unlike the paper [H], so the above correlation 
functions become diagonal. From the horizon crossing relation as 

k = aH= ^e^[/i/2 (3.56) 



we obtain 



c^lnA; = {K'^U + ——]dT 
' 2U dr ' 



2 

m 



'0 

,2 



^ n'UdT. (3.57) 
Then we can calculate the spectral indices as 

d , .0 4 1 fB(2) BU)\ 

_,„«Cf»=--3^(^ + ^j, (3.58) 

« in « SI, »^ 1^ LM4f21±« _ m . (3.59) 



dink K^B{2)\ " "771402(0) + m^02(o) B{2) 

(9 , 5 4 5(4) 

In « 7^,. »= y,, (3.60) 



dink K^B{2Y 
where 

5(277):= 5^ (777^)" 0^(0). (3.61) 

a 

As for the correlation between the adiabatic and the entropic perturbations defined by 

Cab := , ^' ^"' ^ , (3.62) 

we obtain 

^ InC . - - 1^^2^2!!«W±Z^MM (3 63) 

91nfc «:25(0) /€2 5(2) " "7^402(0) +^4^2(0)- ^ 
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Since we succeed in solving the scalar fields evolutions completely, our expressions of the 
spectral indices of the perturbation variables are expressed in terms of the initial fields 
values and the masses of the scalar fields only, unhke the paper |H] where the author did 
not solve the scalar fields evolutions and their expressions of the spectral indices contain 
unknown factors originating from the scalar fields evolutions. 

3.4 Effect of the nonresonant interactions 

In the previous subsection, we considered the free scalar fields. In this subsection, we will 
consider the interacting case. All the interacting systems are classified into the nonresonant 
cases and the resonant cases. By the word "resonance" , we mean all the factors inducing the 
small denominators in the perturbative expansion P]. Therefore the resonance phenomena 
are not confined in the oscillatory dynamical system. In fact, the resonance can also occur 
in the present slow rolling system. In the nonresonant case, the masses of the scalar 
fields do not satisfy any resonant relations as shown in Proposition 7 presented in the 
below. In the resonant case, the masses of the scalar fields satisfy more than one resonant 
relation. In the nonresonant case, the Poincare theorem greatly simplifies the treatment of 
the interactions. In the nonresonant case, all the interaction terms can be removed from 
the evolution equations by performing the suitable transformation of the field variables. 
This process is called the linearization of the evolution equations. Since the linearized, 
that is transformed evolution equations are the free fields equations, they can be solved 
easily. The effects of all the interaction terms are contained in the transformation law of 
the field variables. The fact that the masses of the scalar fields do not satisfy any resonant 
relations guarantees the convergence of the transformation law of the field variables. In 
this way, the evolutions of the interacting scalar fields 0^ can be solved completely in the 
form of the analytic functions of the r function. As for the concrete nonresonant model, 
we determine the evolutions of the scalar fields 0a (t) and the N potential. By using this 
potential, we can calculate the nonlinearity parameters /nl, Qnl representing the non- 
Gaussianity of the Bardeen parameter. These nonlinearity parameters fNLi Qnl are shown 
to be suppressed by the slow rolling parameter. Therefore it is difficult to observe such 
small non-Gaussianity. 

First we present the proposition about the nonresonant condition. In order to express 
the nonresonant condition, we introduce several notations. By a, we represent 

a := (mi,m2, ■ ■ ^m^^) , (3.64) 
k := {ki,k2,- ■ ■,kNs) , (3-65) 

where all ka belong to the nonnegative integer set Zo+: 

ka e Zo+, (3.66) 

Zo+ := {k e Z\k > 0}. (3.67) 

By we represent 

\k\ := (3.68) 

a 

Then the following proposition holds. 



22 



Proposition 7 Suppose that for an arbitrary {a;k) G {1,2,- ■ -jNs} x satisfying 
\k\ > 2, 

{k-a)-aa^O (3.69) 
holds and for an arbitrary a G {1, 2, ■ ■ -, Ns}, 

ml>0 (3.70) 

holds. Then there exists a positive constant 6m satisfying 

\{k ■ a) - aa\ > 5m, (3.71) 

for an arbitrary (a; k) G {1, 2, ■ ■ Ns} x satisfying \k\ > 2. 
For the proof, see Appendix I. 

Under the condition of Proposition 7, we can prove the following proposition. 



Proposition 8 We consider the evolution equation 

^<Pa = -aa(t)a + U(t^), (3.72) 

where fa is the sum of m-th order monomials of (pa (m>2) satisfying 

|/.(0)|</ie^|0r, (3.73) 



where au is the maximum of aa, for < (pQ. Under the condition that a satisfies ^3. 71\ ), 
there exists a positive constant R such that, for (p satisfying \(p\ < R, there exists the 
transformation 

(pa = <^a + Wa{ip), (3.74) 

where Wa is the sum of m-th order monomials of ipa (m>2) and satisfies 

\wam<^W?. (3.75) 

By this transformation law \3. 74\ ), Ij^. 7^ is transformed into the linear differential equation 
as 

d 

— fa = -Oiafa- (3.76) 
OT 

Proposition 8 is well known as the Poincare theorem. [1]] The proof is given in Appendix 
J . This theorem can be applied only to the differential equations containing the linear terms 
plus the small perturbations. Since they are not even the Hamiltonian dynamical systems 
in general, Proposition 8 does not cover the anharmonic oscillator. It is well known that 
the evolution equation of the nonlinear anharmonic oscillator cannot be transformed into 
the linear evolution equation. But for the present purpose of treating the multiple slow 
rolling scalar fields systems, Proposition 8 is enough. All the terms in Wa{'~p) in fl3.74p have 
factors as 

^ (3.77) 



{k ■ a) — aa 
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In order to prove the convergence of this transformation law fl3.74p by the majorant method, 
the nonresonant inequahties fl3.7ip are essential. According to Proposition 8, by (13.741) and 

= 9^a(0) exp [-aaT], (3.78) 

the evolutions of the scalar fields (pa can be expressed as an analytic function of r. 
Next we consider the concrete example whose potential U is given by 

U = ^ml<f)l + ^ml<Pl + ^4>l<Pl, (3.79) 

where ml, m\ are assumed to be nonresonant as in Proposition 7. We obtain 

01 = 0i(O) exp {-mlr) - -£^0i(O)02(o) [exp (-m^r) - exp {-{ml + 2ml)r)}] 

+ (3.80) 

02 = 02(0) exp (-m^r)- ^02(0)0^(0) [exp (-m^r)- exp {-(m2 + 2m2)r)}] 

+ (3.81) 
Then we obtain the potential as 

> = \m + \m - s^^/.mm + ■ ■ ■. (3.82) 

In this model, we consider the nonlinear parameters f^L, g^L [HI] defined by 

5 NabN'^A 

6 {NaN-) 



. 1%^, (3.83) 



where 



We obtain 



25 Nab^N-mN- 
9nl = —; , - - ,3 ' 3.84 

50a(O) 90„(O)a0fe(O) 



72' 



9NL = -^,^^2<^lmm^,-^[^] ■ (3.87) 
144m^ + m2 N'^ U \J^o, 



So we can conclude that the nonlinear parameters /tvl, Qnl are suppressed by the slow 
rolling parameter {m/H^Y. For an arbitrary potential model satisfying the assumptions 
(i), (ii) and (iii) presented above Proposition 6 and for the nonresonant masses as explained 
in Proposition 7, the potential can be written as 



^,N = \Y.^l{Q)+9{m), (3.88) 
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where g is the sum of m-th order monomials of 0a (0) {m > 3) and satisfies 

~gim)<i^c^\m\', (3.89) 

for 10(0)1 < 00- Then we obtain the nonhnear parameters as 

5 1 / m \ / m 



that is, the nonhnear parameters Jml, Qnl are suppressed by the slow rolhng parameters 
3.5 Effect of the resonant interactions 

According to Proposition 8 in the previous subsection, all the nonresonant interaction terms 
can be eliminated by the field transformations 0a — > '■Pa- But in the resonant interaction 
case, such linearization cannot be applied. In this subsection, we calculate the effect of 
the resonant interaction by the iteration method. Perturbatively at least, the resonant 
interactions do not generate any special effects in the N potential. 
As the concrete example, we consider the model defined by 

t^int = A0i0^, ml = 2ml. (3-91) 

It can be solved as 

01 = 0i(O)exp(-m^r) - A0^(O)rexp(-m^r) + • • •, (3.92) 

2A 

02 = 02(0) exp (-mlr) ^0i(O)02(O) Texp {-mlr) - exp {-(m? + mDr)}] 

mf 

+ ■■: (3.93) 

where in the right hand side of 0i the second term proportional to Ar appears because of 
the resonant interaction. The potential is given by 



1 1 ,9.„^ . 1 ,o.„^ 5A 



-N = -02(0) + -02(0) - -—01(0)02(0) + ■ • •. (3.94) 



So the resonant interaction terms generate no singular terms in the A^ potential A^. 



§4 Discussion 

We presented the nonlinear LWL formula and by using it we investigated the evolutionary 
behaviors of the nonlinear cosmological perturbations on superhorizon scales in the universe 
dominated by the multiple slow rolling scalar fields. 

In the excellent study as to the multiple slow rolling scalar fields |5] ] [H ] , the multiple 
scalar fields were decomposed into the adiabatic field and the entropy fields instant by 
instant. By using this decomposition, it was pointed out that the growing of the Bardeen 
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parameter corresponds to the curvature of the trajectory in the scalar fields space. Fur- 
ther in the paper |H ] , the spectral index (the wavenumber dependences) of the Bardeen 
parameter was presented by calculating the differential coefficients at the instant of the 
first horizon crossing. But this study is the local investigation because the author did not 
calculate the S formulae connecting the final amplitudes of the adiabatic and the entropic 
fields variables at the end of the slow rolling phase with the initial amplitudes of them at 
the first horizon crossing, and they calculated only the time derivative of the S formulae at 
the first horizon crossing by using the perturbation evolution equations. So the formulae of 
the spectral indices in the paper ^] have black boxes (unknown factors) which come from 
the S formulae which could not be determined. Therefore as the method for investigating 
the evolutions of the scalar fields perturbations in the long time interval such as in the 
whole slow rolling phase analytically, we present the method using the r function and the 
potential. This method enables us to calculate the S formulae analytically. In fact, our 
expressions of spectral indices in subsection 3.3 can be expressed in terms of the initial 
field values and the masses of the scalar fields only and they have no blackboxes since we 
succeed in calculating the scalar fields evolutions containing the S formulae completely. 
Our method will deepen the understanding of the dynamical evolutions of the multiple 
slow rolling scalar fields in the long time interval with high accuracy. 
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•A Derivation of (12.591) 



In this appendix, we rewrite the momentum constraints ( 12.17^ into the constraints which 
the solution constants of the solution of the locally homogeneous universe must satisfy. 
By using (I2.9p(l2.51i) . the left hand side of the momentum constraint (12.171) is expressed 

as 



1 

a 



H 



a 



(A.l) 



where *f is the Christoffel symbol of 7jj. As for the right hand side of (lA.ip . 



1 diU^-, 



1 a 



'^^'^--2—^ ^^^'^ = -2^^^ 



(A.2) 



and 



%l''dt^ik = -tr 



M-^djMM-^M 



(A.3) 
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and 



and 



Therefore we obtain 



- -D,K 

1 

a 



4a3 



+ - 



a 



dtdi 



to 



a 



Ct + 



1 a 
2^ 



a 



dt—djCT + 2adj 



to 



where we used the solution of M fl2.53l) . 

The right hand side of (12.171) is expressed as 



K 



a 



By summing (12.441) with fluid indices a, we obtain 
[h{pi + Vi)] f = ^y^Cai + ^ I dtaa^ -diPj - -Diuhf - Sadi(t>a 



(A.4) 
(A.5) 



(A.6) 



(A.7) 



+ 0(e3), (A.8) 



where we used Qc 
we obtain 



i)S 



J2a Sadi4>a ^ud as for the scalar fields components 



[ha{(3i + Vai)]s = - XI 't>adi(l)a + 0{€^ 



(A.9) 



We substitute the sum of flA.81) and flA.Qp into h{(3i + Vi) in flA.7l) . 
Then through simple calculations we obtain (I2.59p . 



§B Lie derivative in the leading order of the gradient 
expansion 

The Lie derivatives of the quantity of upper index and the quantity with lower index are 
expressed as 

L{T)X^' = T^'d^X^' -da.T^'X", (B.l) 
L{T)X^ = T"5,X^ + X,9^T^ (B.2) 

respectively. By the Leibniz rule, these definitions are expanded into the tensor of an 
arbitrary rank. For example, the Lie derivative of the metric is given by 

L{T)g,,, = T'^d^g^, + (^^.S^T" + g^^d.T"^. (B.3) 
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From now on, we consider the gradient expansion scheme defined by 

a, = 0(e), r = 0(e), g^o = 0{e), (B.4) 

and the O(e^) order corrections are dropped. For an arbitrary scalar quantity A, in our 
scheme ( 1B.4I) the Lie derivative of A is written by 

L{T)A = T''^A + 0{e^). (B.5) 

By using ( IB. 31) to (700 = — «^ + (3kP^ and (^j-,- = 7ij, we obtain 

L(T)« = T^'^a + a^T^ + 0{e^), (B.6) 
at at 

L(r)7., = T°|7., + 0(e2). (B.7) 
From (IB.Sp and (lB.6p . we obtain 

L{T)h\=T'^h\+0{e% (B.8) 
\ a / at \a I 

From this equation, we can see that the Lie derivative of A/a where A is a scalar quantity 

has the same form as the Lie derivative of the scalar quantity (IB.Sp . By using flB.7l) to 
det(7jj) = a®, we obtain 

L{T)a = T''^a + 0{e^), (B.9) 

therefore we obtain 

LiTn,=T'j^^,, + 0{e'). (B.IO) 

In our gradient expansion scheme, the scale factor a and jij can be regarded as the scalar 
quantity. 



§C Proof of Proposition 1 

The gauge transformation is described by the differential equation (12.701) . So we can 
consider the fi derivative as the gauge transformation. By differentiating (12.701) with respect 
to A, we obtain 

d d^A ^ fd^T\ , ^ fdT\ dA ,,^,d^A , 
M + 2L — — + L(T)-— , (C.2) 



dfi dX^ \ dX^ J \dX J dX dX^ 

and so on. As for D'^A where D is defined by (I2.72p and A is an arbitrary scalar quantity, 
we can prove 

= L(T)D"A, (C.3) 

a/i 
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for an arbitrary natural number n. 
Proof 

If we interpret D^A := A, (1C.3P holds evidently for n = 0. We assume that (1C.3P holds 
for n = k — 1 where k = 1,2, ■ ■ ■. Then for n = k 

—D^A = — S (—-—-— ^D'^'^A 
dji (i/i \ \d\ dX a dt 



d d 
d\ djj, 
'dT_\ 

dx) 



dX \d^ J d dXdfx {d J 

= L{T)D^A (C.4) 

where we use the fact that both D^~^A and the scale factor a have the Lie derivatives of 
the scalar quantity type (IB.Sp . For n = k, (IC.SP holds. By induction, we complete the 
proof. ■ 

By putting A = in ( 1C.3I) and by noticing T(A = 0) = 



■^D"A 
djj, 



= 0. (C.5) 

A=0 



Then we proved that D^-A where A is an arbitrary scalar quantity is gauge invariant. Since 
A/ a where A is an arbitrary scalar quantity has the same Lie derivative as that of the scalar 
quantity A, D^{A/a) is also gauge invariant. In the way similar to the above calculation, 
as for An defined by (12.731) we can prove 



^An = L{T)An. (C.6) 



Then An is gauge invariant. 



§D Proof of Proposition 2 

Since Ci is conserved for arbitrary values of the integration constants C{X = 0), (i with 
C(A = 0) replaced with C(A) is also conserved, that is, Ci is conserved for arbitrary values 
of A: 




(D.l) 



Then by the expression of 6 fl2.80p we obtain 
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Since {d/da)a and {d/d\)a are commutative, we obtain 



d 



d d 



^da J ^ \dadX 
By iterating the same process, we obtain 



Ci 



C =(- 



n-1 



1) 



for n > 2. 



d_d_ 
d\ da 



n Sn 

oa 



0. 



(D.3) 



(D.4) 



§E Proof of Proposition 3 

When P = P{p), 



Then 



a| = -3(, + P). 



'dp\ I fdp_ 
dXj ^1 \da ) „ \ a 
1 1 [dp^ 



3p + P \d\ 
d_ 
dX 



'^Pp + P[p) 



and 



-Ci 

da 



d_d_ 
dXda^ ^ 
d_ 
dX 



sJ'^Pp + Pip) 



0. 



(E.l) 



(E.2) 



(E.3) 



since {d/dX)a and {d/da)a are commutative. Then we can use Proposition 2, we complete 
the proof. 

§F Proof of Proposition 4 



For = 0, by solving (13. lip f l3.12p we obtain 

|p_e-3A^p(0)| < 

10-0(0)1 < v^N. 



[F.l) 

:f.2) 



If p{0), 0(0) ~ 1, Proposition 4 is right for A; = 0. We assume that Proposition 4 holds for 



k = 0,1,2,- ■ ■,k-l: 



( ^ 



dX^ 



p 



SI 



(^ = 0,1,2,---,A;-1), 



(F.3) 
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therefore 



if 



|-) mp,CT,N) 



61, (2 = 0,l,2,---,fc-l), 

\f{c^,p,CT,N)\<l, 



(F.4) 
(F.5) 



for a proper complex domain containing the real interval which we consider. For k, the 
evolution equations are 







-1 




± 


/ Qk 


dN 





p 



V 

-3 



Qk 



p + r] 



Qk 



2xk 



Qk 



p + 1] 5_ 



Qk 



(F.6) 

p + v'St (F.7) 



where all the coefficients bounded by positive constants are omitted except —3 in ( IF. 71) . 
By solving (iradRTl) . for < < l/r/^, we obtain 



Qk 



Qk 



< 



Qk 



9A^- 



0(0) 



Qk 



P 



< e 



-3N 



d' 



dX 



kP(^) 



dX' 



+ r] 



:P{0) 



Qk 



dX^ 



m 



V 



Qk 



dX' 



:P(0) 



+ 6' 



(F.8) 
(F.9) 



Therefore for the initial conditions ( ]3.17p . for k, (]3.18p(]3.19p hold. By induction, for all 
nonnegative integers k, Proposition 4 holds. 



§G Proof of Proposition 5 



First we consider A; = case. By the mean value theorem, we obtain 

d 



where 



(G.l) 
(G.2) 



\f\<\p\+v'<e-''' + V^ 
where we used Proposition 4 in the last inequality. By solving (IG.ip . we obtain 

I A0I < exp {r]'^N) (| A0(O) 1+7]"^ + r]^N) . (G.3) 

Therefore for < N < l/r]"^, under the initial condition fl3.22p . (13.211) holds for k = 0. 
Next we consider a positive integer k case. In this case, we notice the following fact. By 
Proposition 4, if an complex analytic function / satisfies 



\f{<P,p,CT,N)\<l, 



then 



Qk 



f{(f),p,CT,N) 
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<6l 



(G.4) 
(G.5) 



By differentiating (1G.1|) with respect to A and by using the above fact and Proposition 4, 
we obtain 

^ ^ ^ A<P = V'(^] A0 + v'S'^ + v') , (G.6) 



dN V ^A'^ 



which is solved as 



A0 



< exp {r{^N) 



gk 



9A^' 



:A0(O) 



(G.7) 



Therefore for < < l/?7^, under the initial condition fl3.22p . (13.211) holds for a positive 
integer k. 



•H Proof of Proposition 6 



Lemma 1 



ba{0) exp [-m^r]+ < /ic^oexp [-m^r] > 



Proof We consider the evolution equation 

^(f)a = -ml(j)a + /a(0), 

where /a(0) is the sum of m-th order monomials (m > 2) satisfying 

|/a(0)|<mV^l0r, 



for 1 01 < 00- From (lH.2p . we obtain 

d_ 

Ft 

which is solved as 



^m<-m'\6\ + ^\6\\ 



< 00 exp [-mV](l + 0(/i,)). 



(H.l 

(H.2 
(H.3 

(H.4 
(H.5 

(H.6 
(H.7 

(H.8 
(H.9 

I /a (0)1 < mVc0oexp [-mV], (H.IO 

where we used (IH.Sp . By solving the above evolution equation under the above initial 
condition, we obtain 

|A0| < O(/ic)0oexp [-mV]. (H.ll 

We complete the proof. ■ 



We consider 0a satisfying 



—0a = -mlcpa, 



The difference A0a : = 



d_ 

0a(O) = 0a(O). 

ha satisfies 



9 ^ . 

A0a(O) = 0, 



-m^A0a + /a(0), 



where 
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Lemma 2 



9' 



Qk 



<^ I v^a - ^v>a(0) • exp [-m^r]+ < ^^/i^^oexp [-mV] > . (H.12) 



Proof As for the evolution equation as 



d f d' 



Or \dX 



Qk 



where 



Qk 



-k] 



< m 



2 /^c 



Qk 



4>o 



Qk 



(H.13) 



2 A'-c c-fci i |2 



00 



< m He exp [— m r] 



+ 6^iJ,cm^(f)oexp [-2mV].(H.14) 



We perform the calculations similar to the proof of Lemma 1. We complete the proof. 



Lemma 3 



f Qk\ ( ^ \ ( ^ \ ^ 



(H.15) 



Proof We use the Leibniz rule for the A derivative and use Lemma 1, Lemma 2. We 
complete the proof. ■ 



Lemma 4 For n> 1 



A{2n,Q) = Y,{KT€ + F2n{ci>) 



where F2n{4>) is the sum of m-th order monomials of (p (m>3) and satisfies 



|i^2n(0)| </icm^"-|<^l^ 
00 



for 101 < 00- 



(H.16) 



[H.17) 



Proof We use (lH.2p . We complete the proof. 
As for 

N = K^ dr 
Jo 

and 



Qk \ .00 / Qk 



9A^ 



kl 



and 



A{0,k) = 4 / dr 



Qk 



Qk 



(H.18) 
(H.19) 
(H.20) 



for n > 1 where A(2n,0) is expressed as a polynomial of 0(r) ( lH.16p . we use Lemma 3. 
Then we complete the proof of Proposition 6. 
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§1 Proof of Proposition 7 

By a^, we mean 



Then we obtain 
Therefore for k satisfying 



Om '■= niax{m^}, (I.l) 

a 

am := min{m^}. (1.2) 

a 

\{k ■ a) - aa\ >\k\am- aM- (1-3) 



\k\ > ^i±^ =: L, (1.4) 



where 6i is a positive constant, we obtain 

\ik ■ a) - aa\ > 6i. (1.5) 

Since the number of k satisfying \k\ < L is finite, there exists a positive constant 62 
satisfying 

\{k-a) - aa\ > S2, (1.6) 
for \k\ < L. Therefore for 5m defined by 

6m ■■= mm{6i,62}, (1.7) 

the inequahty fl3.7ip holds. We complete the proof. 

§ J Proof of Proposition 8 

By replacing (pa with 0o0a and by replacing r with t/ ficC^M, we can assume that |/a(0) | < 1 
for 101 < 1. In order that the transformation law fl3.74p can reduce the original equation 
fl3.72p to the exactly linear equation fl3.76p . Wa must satisfy 

dw ~ 

- —-^abfb + aaWa = faif + w). (J.l) 

We decompose Wa into the sum of the i^-th order polynomials: 

wa = j2 ^a, < = E ^aV^ ( J.2) 

K>2 \k\=K 

where 

if'' := ifl'--- cp'^l', \k\:=k^ + --- + k^s- (J.3) 
We decompose ( IJ.ip into i^T-th order part as 

c)ni^ ~ 

\k\=K 
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for K = 2, ■ ■ -. We mean the coefficient of if'' in by Then c?^ (|fc| = K) is the 
polynomial of < K) and < K — 1) with integral coefficients and can be 

expressed as 

-{k ■ a) + aa 

From now on, we will prove the convergence of the transformation Wa by the method of 
the majorant. F{ip) defined by 

oo 

F{^) = Y,{^, + --- + ^^^f (J.6) 

K=2 

is the majorant of fa- We consider the equation as 

6m-W{v) = F{^ + W{^)), (J.7) 

where Sm is defined in Proposition 7, Since F{f) is a independent, Wa{(p) is a independent. 
Then Wa{(p) is simply written as W{(f). We will prove that W{ip) is the majorant of Waii^). 
By expanding W{ip), F{(p) F{ip + W(ip)) as 

W{^) = W^V, F{¥^) = E ^V, F{V + W{^)) = E DV', (J-8) 

|fc|>2 |fc|>2 |A;|>2 

we obatin W'' = /5m- For \k\ = 2, = F^ , so \w^\ < W^. We assume that \w^\ < 
for \k\ < K — 1. Since c?^ {\k\ = K) is the polynomial of (|/c| < K) and (|A;| < 
K — 1) with integral coefficients and is given by substitution f^^F^^w^^ in 
the polynomial representing d^, < {\k\ = K). Therefore we obtain < W'' 
{\k\ = K). By induction, W{ip) is the majorant of Wa{'~p)- As for the right hand side of 

(IlrD, 



so long as 1$ + NsW{ip)\ < 1. So (]J.7p can be written as 

iVs(iV5 + 5m)W^ - [5m - {2Ns + 6^)^^ + = 0, (J.IO) 
whose solution is given by 
1 



W 



2NsiNs + 5m) 



6m - {2Ns + 5m)^ - ^5l-25m{2Ns + Sm)^ + Sl^^ (J.ll 



for 1$ + NsW{(p)\ < 1, using the fact that W = for ipa = 0. Since W is the majorant of 
Wa, we can prove the convergence of Wa- 

This proof in Appendix J is based on [21] ] . 
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